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Abstract. It is almost 80 years since Hans Bethe described the level density as a non-interacting gas 
of protons and neutrons. In all these years, experimental data were interpreted within this picture of 
a fermionic gas. However, the renewed interest of measuring level density using various techniques calls 
for a revision of this description. In particular, the wealth of nuclear level densities measured with the 
Oslo method favors the constant-temperature level density over the Fermi-gas picture. From the basis of 
experimental data, we demonstrate that nuclei exhibit a constant-temperature level density behavior for 
all mass regions and at least up to the neutron threshold. 


PACS. PACS-key 21.10.Ma, 25.20.Lj, 25.40.Hs 

1 Introduction 

The level density of atomic nuclei provides important in¬ 
formation on the heated nucleonic many-particle system. 
In the pioneering work of Hans Bethe in 1936 pQ, the level 
density was described as a gas of non-interacting fermions 
moving in equally spaced single-particle orbitals. Despite 
the scarce nuclear structure knowledge at that time, this 
Fermi-gas model contained the essential components apart 
from the influence of the pairing force between nucleons in 
time-reversed orbitals. The understanding of such Cooper 
pairs and the nuclear structure consequences, was first 
realized twenty years later through the Bardeen-Cooper- 
Schrieffer (BCS) theory [ 2 ]. 

Unfortunately, the new knowledge of the pairing force 
was implemented in an oversimplified way by still assum¬ 
ing a gas of fermions. The original Fermi-gas level density 
curve was simply shifted up in energy by A and 2 Z\ for 
the description of odd-mass and even-even mass nuclei, 
respectively. Later, this energy correction was found to be 
too large, and the shift was somewhat back-shifted again. 
The resulting back-shifted Fermi-gas model has since then 
been very popular and has been the most common descrip¬ 
tion of nuclear level densities for decades [3|. These ma- 
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nipulations maintained the typical exp (2 y/aE) Fermi-gas 
level density shape with excitation energy E. 

There is a growing interest in the nuclear science com¬ 
munity for the study of nuclear level densities. The in¬ 
troduction of novel theoretical approaches and fast com¬ 
puters has opened the way for a microscopical descrip¬ 
tion of heavier nuclei up to high excitation energy |4 . 
Thus, experimental level densities represent a basic testing 
ground for many-particle nuclear structure models. They 
also have an increasing importance for various nuclear ap¬ 
plications. Nuclear level densities play an essential role in 
the calculation of reaction cross sections applied to astro- 
physical nucleosynthesis, nuclear energy production and 
transmutation of nuclear waste. 

The number of levels is exponentially increasing with 
excitation energy and typically for the rare earth region, 
the level density increases by a factor of one million when 
going from the ground state up to the neutron binding 
energy. This tremendous change has implications for the 
theoretical interpretations since levels start to overlap in 
energy at high excitation energies, typically around 10 - 
15 MeV. Thus, nuclear excitations are often divided into 
three regimes [5], classified by the average level spacings 
D compared to the 7 -decay width T. These key numbers 
are connected to the level density and life time by p — 
1/D and r = %/T, respectively. Figure [l] illustrates the 
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Fig. 1 . (Color online) The energy regimes of nuclear excita¬ 
tions. 


connections between the three regimes i.e. the discrete, 
quasi-continuum and continuum regions. 

Experimentally, the quasi-continuum region is often 
defined as the region where the experimental resolution 
AE res prevents the separation into individual levels. This 
appears when AE res > D. One should be aware that this 
definition depends on the selectivity of the reactions used 
and the specific experimental detectors and conditions. In 
this work, the experimental quasi-continuum region starts 
when AE res « D « 100 keV, i.e. when the level density 
exceeds p & 10 MeV -1 . 

Today various experimental techniques are used to de¬ 
termine nuclear level densities. The obvious method in the 
discrete region of Fig. [l] is to count the number of levels 
AN per excitation energy bin AE , readily giving the level 
density as p(E) = AN/AE. This requires that all levels 
are known. Making a short review on known discrete lev¬ 
els, e.g. in the compilation of Ref. 0, we find that the 
data set starts to be incomplete where the first Cooper 
pairs are broken, i.e. at excitation energies E « 2 A. This 
is simply because the high density of levels makes con¬ 
ventional spectroscopy difficult. A rule of thumb is that 
when p & 50 — 100 MeV -1 , the number of missing levels 
is dramatically increasing. 

Level density information at the neutron separation 
energy S n can be extracted from neutron capture reso¬ 
nance spacings D [7]. The spin selection of these reso¬ 
nances depends on the target ground-state spin and the 
neutron spin transfer, giving level densities for a narrow 
spin window, only. In particular, thermal neutrons with 
angular momentum transfer i = 0 results in a spin window 
of ^target ±1/2. At higher excitation energies the method 
of Ericson fluctuations can be used |5 . A well known tech¬ 
nique |8| is to extract the level density from particle evap¬ 
oration spectra. The measurement must be carried out 


at backwards center-of-mass angles to avoid direct reac¬ 
tion contributions. The method also requires well deter¬ 
mined optical potential model parameters. Furthermore, 
the level density can be extracted in the two-step cas¬ 
cade method [9] from simulations provided that the q-ray 
strength function is known. An upcoming technique, see 
e.g. Ref. [TO] , is to use high-energy light-ion reactions like 
(p,_p 7 ), (e, e') and ( 3 He, t) at small angles with respect to 
the beam direction. These reactions select the population 
of discrete levels with certain spin/parity assignments, in¬ 
formation that is essential in the detailed understanding 
of level densities. This method requires that the spacings 
of the selected levels are larger or comparable to the ex¬ 
perimental particle resolution. 

In the present work we will show level densities of sev¬ 
eral nuclei measured at the Oslo Cyclotron Laboratory 
(OCL). The technique used, the Oslo method |111112] . is 
unique in the sense that it allows for a simultaneous de¬ 
termination of the level density and the q-ray strength 
function without assuming any models for these functions. 

The structure of the manuscript is as follows: Section 
2 describes briefly the experimental set-up and the Oslo 
method, and previously measured level densities will be 
discussed in Sect. 3. Finally, concluding remarks and out¬ 
looks are presented in Sect. 4. 


2 The Oslo Method 

The Oslo method is based on a set of q spectra as a func¬ 
tion of excitation energy. The excitation of the nucleus is 
performed by light ion reactions, e.g. (d,pq), (p,p'q) and 
( 3 He, cry) where the energy of the charged ejectile deter¬ 
mines the excitation energy. 

A schematic drawing of the set-up is shown in Fig. [2] 
A silicon particle detection system (SiRi) [13], which con¬ 
sists of 64 telescopes, is used for the selection of a certain 
ejectile type and to determine its energies. The front AE 
and back E detectors have thicknesses of 130 /im and 1550 
/im, respectively. SiRi is usually placed in backward angles 
covering 0 = 126° to 140° relative to the beam axis. Co¬ 
incidences with q rays are performed with the CACTUS 
array consisting of 26 collimated 5” x 5” Nal(Tl) de¬ 
tectors with a total efficiency of 14.1% at E 1 = 1.33 MeV. 

Figure [3] demonstrates the previously measured 
237 Np(d,pq) 238 Np reaction [15] and how one can proceed 
from the raw particle-q coincidences to the first genera¬ 
tion or primary q-ray spectra. The first step is to sort the 
events into a raw particle-q matrix R(E, E 7 ) with proper 
subtraction of random coincidences. Then, for all initial 
excitation energies E, the q spectra are unfolded with the 
Nal response functions giving the matrix U{E^E 1 ) 16 . 
The procedure is iterative and stops when the folding T 
of the unfolded matrix equals the raw matrix within the 
statistical fluctuations, i.e. when J-(U) ~ R. 

The primary q-ray spectra can be extracted from the 
unfolded total q spectra U of Fig. [ 3 ] (b). The primary q 
spectrum at an initial excitation energy E is obtained by 
subtracting a weighted sum of U(E'^E 1 ) spectra below 
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Fig. 2. (Color online) Typical particle -7 coincidence set-up for 
the Oslo method. The 64 silicon particle telescopes of SiRi are 
placed in the vacuum chamber at the center of CACTUS. 


excitation energy E : 

P(E, E y ) = U(E, E y ) - Y, W ( E ’ E ')U{EE 7 ). (1) 

E'<E 

The weighting coefficients W(E,E f ) are determined in an 
iterative way described in Ref. m- After a few iterations, 
W(E,E f ) converges to P(E, E y ), where we have normal¬ 
ized each 7 spectrum by ^ E P(E,Ey) = 1. This equal¬ 
ity P « W is exactly what is expected, namely that the 
primary 7 -ray spectrum equals to the weighting function. 
The validity of the procedure rests on the assumption that 
the 7 -energy distribution is the same whether the levels 
were populated directly by the nuclear reaction or by 7 de¬ 
cay from higher-lying states. This is illustrated in Fig. [4j 

The statistical part of this landscape of probability, 
P(E , £ , 7 ), is then assumed to be described by the product 
of two vectors 

P(E,E>p(E-E 7 )r(^), (2) 

where the decay probability should be proportional to the 
level density at the final energy p(E — E 7 ) according to 
Fermi’s golden rule jT9[f20] . The decay is also proportional 
to the 7 -ray transmission coefficient T, which according 
to the Brink hypothesis [21], is independent of excitation 
energy; only the transitional energy E 1 plays a role. 

The values of the vector elements of p and T are varied 
freely to obtain a least-% 2 fit to the experimental P ma¬ 
trix. Figure [ 5 ] shows that the product of p{E — E 7 )T(E 7 ) 
of Eq. Q fits very well the primary spectra of 238 Np at the 
six excitation energy bins shown. There are actually four 
times more spectra and the same quality of fit is obtained 
for all spectra. 

However, the nice fit to the data does not mean that p 
and T are uniquely determined. There are infinitely many 
functions that make identical fits to the experimental P 



E y (MeV) 


Fig. 3. (Color online) Initial excitation energy E versus 7 - 
ray energy E 1 from the 237 Np(d,p 7 ) 238 Np reaction [15]. The 
raw 7 -ray spectra (a) are first unfolded by the Nal response 
function (b) and then the primary 7 -ray spectra are extracted 
(c). 


matrix. These p and T functions can be generated by the 
transformations m 

p(Ei - E 7 ) = Aexp[a(Ei - E^)\ p(E { - E . 7 ), (3) 

T(E^) = B exp(aE^)T(Ej). (4) 

Thus, the normalization of the two functions requires in¬ 
formation to fix the parameters A , <a, and H, which is not 
available from our experiment. 
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Fig. 5. (Color online) Primary 7 -ray spectra from the 237 Np(d,p 7 ) 238 Np reaction [15] gated on various initial excitation 
energies E (crosses). The spectra are compared to the product of the level density and transmission coefficient functions 
i.e. p(E — E 7 )T(E 7 ) (blue lines). The excitation energy gates are 121 keV broad. 





The same 


y-energy 

distribution 



Fig. 4. (Color online) The assumption of the extraction of pri¬ 
mary 7 spectra is that we obtain the same energy distribution 
(yellow) whether the levels of the energy bin E 2 are populated 
directly by e.g. the (d, p) reaction or by 7 transitions from ex¬ 
citation bin Ei. 


Fig. 6. Normalization procedure for the experimental level 
density. The experimental data are normalized to known dis¬ 
crete levels at low excitation energy and to the level density 
extracted at S n from neutron capture resonance spacings Do. 
The two set of arrows indicate where the data are normalized 
by means of Eq. <§• 


For the level density, we use two anchor points to de¬ 
termine A and a. The procedure is demonstrated for the 
case of 238 Np in Fig. [6] The level density is normalized 
to known discrete levels at low excitation energy. At high 
excitation energy, we exploit the neutron resonance spac¬ 
ings Dq at the neutron separation energy S n . In order 
to translate the spacing to level density, we use the spin 


distribution El 

g(E, I) ~ 2 ^J~ 2 1 exp [-(/ + l/2) 2 /2 a 2 ] , (5) 

where E is excitation energy and I is spin. Here, the spin 
cut-off parameter a is rather uncertain, and may intro- 
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Fig. 7. (Color online) Experimental level densities obtained with the Oslo method. For references, see 


Table [I] 


duce an uncertainty of a factor of 2 in p(S n ). For cases 
where Dq is unknown, one has to use systematics to esti¬ 
mate p(S n ). Since the 7 strength is not discussed in this 
work, the determination of the B parameter is irrelevant. 
Further description and tests of the Oslo method are given 

in Refs. inma. 


3 Experimental level densities 

The Oslo group has now published more than 70 level den¬ 
sities. Some of these are for one and the same nucleus, but 
populated with different reactions. The full compilation of 
level densities (and 7 -ray strength functions) is found in 
Ref. [3C|. We will here focus on some selected nuclei in 
order to convey to the reader a general picture of how the 
level density behaves as function of mass region and exci¬ 
tation energy. The references and key numbers are given 
in Table |T| 

There is a tremendous span in level densities when 
going from mass around A « 40 to 240. Figure [7] shows 
that the level density at excitation energy E « 6 MeV for 
238 Np is almost ten million times higher than for 46 Ti. It 
seems obvious that more nucleons A produce higher level 
density, but this conclusion is not strictly true. 

Indeed, we see that the level density for 90 Y is larger 
than for 89 Y. Also 117 Sn has more levels than 116 Sn. How¬ 
ever, 45 Ti has more levels than 46 Ti and the same is true 
for 163 Dy versus 164 Dy. The pieces fall into place if pair¬ 
ing is taken into account: For neighboring isotopes, the 
nucleus with either one or more unpaired protons or neu¬ 
trons displays 5 — 7 times more levels. An increase in A 


only gives negligible effects between neighboring nuclei, as 
e.g. found for the even-even 160 462 ,i 64 j}y i so t 0 pes [51] . 

The number of quasi-particles is the basis building 
blocks for the level density. However, the number of active 
particles making the level density is not that different for 
238 Np and 46 Ti. The odd-odd 238 Np has about 6 quasi¬ 
particles at E « 6 MeV, while 46 Ti has about 2 quasi¬ 
particles. The second factor needed to understand the gen¬ 
eral behaviour of the level density is the density of avail¬ 
able single-particle orbitals in the vicinity of the Fermi 
level where the quasi-particles can play around. In sum¬ 
mary, the large span in level density between A = 40 and 
240 is a combined result of the number of quasi-particles 
and the density of near lying single-particle orbitals. 

Another striking feature of Fig. [7] is that, above a cer¬ 
tain excitation energy E , the level densities are linear in 
a log-plot, which corresponds to an exponential function. 
The 45 Ti has an exponential behavior for E > 4 MeV, 
whereas 238 Np has it for E > 1 MeV. Generally, the ex¬ 
ponential behavior sets in when the first nucleon pairs 
are broken, i.e. E > 2A. With a pairing gap of A « 
12 A -1 / 2 MeV, we find for the two cases E = 3.6 and 
1.6 MeV, respectively. The level density is not only de¬ 
pendent on the atomic mass number A, but also on the 
number of unpaired nucleons around the ground state. 
As an example, the exponential behavior is delayed for 
even-even nuclei compared to odd-even nuclei. This is e.g. 
apparent for 116,117 Sn and 163 d64j7y 

The interpretation of the pure exponential growth above 
E > 2A has recently been discussed in Ref. m- For the 
micro-canonical ensemble, the entropy is related to the 
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Table 1. Parameters to describe the level density. 


Reaction 

Sn 

(MeV) 

Tct 

(MeV) 

Eg 

(MeV) 

p(Sn) 

(10 6 MeV 1 ) 

Refs. 

(p, P r° Ti 

9.530 

1.55 

-2.33 

0.0014(7) 

m 

(d,p) 46 Ti 

13.189 

1.70 

-2.07 

0.0047(10) 

m 

(p,pT y 

11.478 

1.00 

0.48 

0.060(20) 

m 

(d,pf° Y 

6.857 

1.00 

-1.38 

0.0038(8) 

m 

( 3 He,a) 116 Sn 

9.563 

0.76 

-0.03 

0.40(20) 

25,26,27,28 

( 3 He, 3 He') 117 Sn 

6.944 

0.67 

-0.40 

0.091(27) 

25,26,27,28 

( 3 He,a) 163 Dy 

6.271 

0.59 

-1.55 

0.96(12) 

m 

( 3 He, 3 He') 164 Dy 

7.658 

0.60 

-0.66 

1.74(21) 

m 

(d,p) 238 Np 

5.488 

0.41 

-1.35 

43.0(78) 

m 


level density by 


S(E) * In p(E), 

(6) 

with 


T{E) = [dS(E)/dE]~ 1 . 

(7) 


Thus, the temperature is constant when In p(E) is linear in 
E. The common expression for the constant-temperature 
level density is [5] 


Pct{E) = —— exp 
J- CT 


E — Eq 
Tct 


( 8 ) 


where Tqt is determined by the slope of In p(E). This 
is the key characteristic of a first-order phase transition. 
Our interpretation of the data is that the energy goes into 
breaking Cooper pairs and thus the temperature remains 
constant. The obvious analogy is the melting of ice to 
water at constant temperature. The fit parameters Tqt 
and Eq for the nuclei shown in Fig.[7]are listed in Table[l] 

Figure [8] shows that the constant-temperature level 
density model fits the data of 164 Dy much better than the 
Fermi gas model. The reduced x 2 fit value evaluated be¬ 
tween excitation energies E = 1.9 and 6.2 MeV is about 
a factor of 10 better for the CT model. One should of 
course mention that the CT model is not a perfect model. 
There are some oscillations in the experimental log-plot 
mainly due to the onset of the two and four quasi-particle 
regimes. 

Recommendations are available for the Fermi-gas spin- 
cut off parameter a (A, E) [32]. However, for the CT-model 
there are no clear recommendations. According to the 
original expression of Ericson [5]: 


a 


2 


OT 

~w 


(9) 


we see that for the CT model, only the moment of in¬ 
ertia O will depend on E, since T = Tct is a constant. 
There are indications |33] that the rigid moment of iner¬ 
tia is reached at the neutron separation energy, i.e. 0(E = 
Sn) = ©rigid- However, how to parameterize a to lower ex¬ 
citation energies is still uncertain. In addition to the uncer¬ 
tain spin distribution, the parity distribution is important 
for lighter nuclei and nuclei close to magic numbers. There 
are only few experimental data available in order to pin 


down the spin and parity distributions and this represents 
perhaps the most important issue for future research on 
nuclear level densities. 

The pair breaking process will continue at least up 
to the binding energy, but probably to much higher exci¬ 
tations. In several evaporation experiments the constant 
temperature behavior has been observed up to E « 15 MeV 
[341135] . However, at the higher excitation energies the 
temperature will start increasing with excitation energies, 
and the Fermi gas description m will to an increasing 
degree describe the data. 


4 Concluding remarks and outlook 

The present work describes the Oslo method and results 
using light-ion reactions on stable targets. The method 
is demonstrated by means of the 237 Np(d,p 7 ) 238 Np reac¬ 
tion, which provides almost ideal conditions for the Oslo 
method. One important aim of these investigations is to 



Fig. 8. (Color online) Experimental level densities of 164 Dy 
compared with the constant-temperature (CT) and the Fermi 
gas (FG) level density formulas. The reduced least-y 2 fit is best 
for the CT model. 
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obtain more predictive power for extrapolation of the level 
densities to nuclei outside the valley of /^-stability. 

From light ( 45 Ti) to the heavy elements ( 238 Np) we 
find certain common features of the level densities. It turns 
out that the pairing force is responsible for the mecha¬ 
nism behind these features. The level density above the 
excitation energy of 2 A behaves as an exponential curve 
characterized by the nuclear temperature. Furthermore, 
neighboring nuclei have the same temperature and their 
respective level densities appear parallel in a log plot. 
Odd-mass nuclei exhibit about 5 — 7 times more levels 
than their even-even partners. There are indications that 
each unpaired nucleon carry approximately an entropy 
of 1.6 — 1.9corresponding to the multiplicative factor 
above. 

In the field of experimental level densities, there are 
several skeletons in the closet. In particular, it is crucial 
to get more information on the true spin and parity dis¬ 
tributions as function of excitation energy and mass. The 
few experimental data on these quantities also have im¬ 
plications for applications when e.g. calculating reaction 
rates. Today, such information is often taken from theo¬ 
retical models without verifications based on experimental 
observations. 

In the energy regime from the ground state up to the 
energy where the breaking of the first Cooper pairs set in, 
the level density depends strongly on the specific nuclear 
structure properties of the nucleus studied. For even-even 
nuclei, the level density is determined by the collective 
modes like rotation and vibrations. For the odd-mass nu¬ 
clei the unpaired valence nucleon already carries a quasi¬ 
particle entropy that adds to these collective degrees of 
freedom. 

A great challence is to obtain data also for nuclei far 
away from the /^-stability line, nuclei that are essential for 
nucleo-synthesis networks and other applications. Since 
it is a rather uncertain approach to extrapolate results 
from more stable nuclei to these nuclei, novel experimen¬ 
tal techniques are needed. The Oslo group has already 
initiated experiments to measure level densities for short¬ 
lived nuclei. The so-called /3- Oslo method is based on the 
population of rare nuclei by means of /3 decay. The first 
experiment, which was performed at the National Super¬ 
conducting Cyclotron Laboratory at Michigan State Uni¬ 
versity (MSU) using the SuN detector has already been 
published [36] . and several new experiments have been 
carried out. Also inverse kinematic experiments with ra¬ 
dioactive beams are planned at HIE-ISOLDE. 

The constant-temperature level density behavior is ex¬ 
pected also for nuclei far from the /^-stability line. How¬ 
ever, when the particle separation energies approaches the 
paring energy of E = 2A, the individual low-lying discrete 
levels will play an increasing role. Furthest away from the 
valley of stability, discrete spectroscopy of low-lying levels 
will probably be the way to proceed. In such scenarios the 
concept of level density will lose its importance and be 
replaced by discrete spectroscopy. 
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